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Abstract 

We report on an experiment to detect non-classical correlations in a highly mixed state. The 
correlations are characterized by the quantum discord and are observed using four qubits in a liquid 
state nuclear magnetic resonance quantum information processor. The state analyzed is the output 
of a DQC1 computation, whose input is a single quantum bit accompanied by n maximally mixed 
qubits. This model of computation outperforms the best known classical algorithms, and although 
it contains vanishing entanglement it is known to have quantum correlations characterized by the 
quantum discord. This experiment detects non-vanishing quantum discord, ensuring the existence 
of non-classical correlations as measured by the quantum discord. 



PACS numbers: 03.67.Lx, 76.60.-k. 



Entanglement has often been synonymous with quantum. Over the past several decades, 
entanglement has been well studied in the context of quantum information and computation 
as it plays a crucial role in many protocols, such as teleportation and superdense coding pQ . 
It has been shown as a vital component to the speedup exhibited by pure state quantum 
computation |2J. However, in the context of mixed state quantum computation it is not 
known what role entanglement plays. In fact, there is a computational model that con- 
tains limited entanglement [3J, yet an advantage over current classical methods appears to 
exist |4J. This model of mixed state quantum computation is known as DQC1, or determin- 
istic quantum computation with one quantum bit, and has been of great interest in recent 
years. 

While DQC1 contains very little, or no entanglement, it does contain non-classical cor- 
relations as measured by the quantum discord [5]. Quantum discord is a measure of the 
correlations that exist in excess of those present in classical states. It is measured by the 
difference in two classically equivalent formulations of the mutual information, where a non- 
zero value indicates a deviation from purely classical correlations [SI [7j. It is not yet known 
whether or not quantum discord assists quantum algorithms, but it is a good candidate for 
the computational advantage offered by the DQC1 model, and understanding it better will 
undoubtably provide insights into the workings of quantum systems and algorithms. While 
it has been shown that almost every quantum state has quantum discord [TU], Datta et al. 
showed that on average, the quantum discord present in a DQC1 algorithm drops with a 
decrease in polarization [5J. Typical liquid state nuclear magnetic resonance (NMR) exper- 
iments are performed at room temperature, and the unentangled initial states have very 
small polarization. Therefore, the question remains, is it possible to experimentally detect 
quantum discord in a DQC1 algorithm where the polarization is very small? 

We report on the experimental detection of non-classical correlations in the output state 
of a DQC1 algorithm in NMR using a state-independent non-zero discord witness. Similar to 
entanglement, calculation of the quantum discord requires full tomographic data, therefore 
witnesses that detect the existence of these correlations are a practical alternative. Before 
describing our results in detail, we define both the DQC1 model and quantum discord. We 
then explain the method used to detect non-classical correlations, followed by a description 
of the physical system. Finally, we conclude by discussing the significance of these results. 

DQC1: Deterministic quantum computation with one quantum bit (often misnamed 
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one clean qubit) is a model of computation where only one qubit is slightly polarized away 
from the maximally mixed state. The DQC1 circuit seen in Figure [T] calculates the real and 
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FIG. 1. The DQC1 circuit where the top qubit has a bias of e towards the ground state. Mea- 
surements of (cr x ) and {a y ) yield the real and imaginary parts of eTr(C/ n )/2 n , respectively. In the 
initial state, the identity terms are recognized as properly normalized states such that tr(p) = 1. 



imaginary components of the trace of the unitary U n upon measurement of the expectation 
values of Pauli operators a x and a y , respectively. Estimating the trace of a unitary matrix 
is a problem with no known efficient classical algorithm, thus, a DQC1 quantum computer 
outperforms its current classical counterparts. Interestingly, there is never entanglement 
between the top qubit and the rest of the system. With any other bipartite split there 
can be a small amount of entanglement that does not grow with the size of the system [3] , 
suggesting that entanglement may not be the cause of the quantum speedup in DQC1. 
Despite the small amount of entanglement found in this model, it contains a significant 
amount of non-classical correlations as measured by quantum discord for a randomly chosen 
unitary [5]. 

Discord: Quantum discord measures non-classical correlations in bipartite states that 
are not necessarily captured by entanglement [6j[7]. These non-classical correlations include, 
but are not limited to, entanglement, and are measured by the difference of two expressions 
for the mutual information, I (A : B) := H(A) + H{B) - H(A, B) and J(A : B) := H{B) - 
H(B\A). For states described by a classical probability {pi}, H(A) = —J2iPi^°SPi * s the 
Shannon entropy, and H(B\A) is the conditional entropy, averaging the entropy of system 
B over the different outcomes on system A [8]. If A and B refer to quantum mechanical 
systems, then their states are described by operators on Hilbert spaces Hab = Ha ® H& 
and H(A) = — Tr(pyi log p^) is the von Neumann entropy. Classically the two expressions 
for the mutual information are equivalent, I(A:B) = J(A:B). However, if all possible 
measurements on a subsystem of the bipartite state disturb the quantum state in Hab, then 
there exist correlations that are not found in completely classical systems, and I {A : B) ^ 
J (A : B). The quantum discord D(A : B) is defined as the minimum difference between the 
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two formulations of the mutual information, and reduces to 

D(A : B) = H(p A ) - H(p) + mm {Ek} ^p fc #(p B | fc ) , (1) 

k 

where pa = Tr#(p) is the reduced density matrix of system A, {E k } is a complete set 
of orthonormal projectors on "H^ such that Ylk — \ Pk is the probability of observing 
outcome k on system A, and ps\k ■= T^A((Ek <8> Ib)p) is the state of system B conditional 
on the measurement of system A returning measurement outcome k. 

Employing Equation [T], the amount of discord in a given state can be estimated relying on 
a full description of the bipartite state. In a recent experiment, full state tomography from a 
two qubit optical implementation of DQC1 was used to estimate the quantum discord [U] . On 
the other hand, to detect non-zero quantum discord, there exist methods that do not require 
full state tomography The experiment presented in this letter follows one such proposal of 
a state independent non-zero discord witness [IT] . The procedure shows that by writing the 
state of the bipartite system as pab = Ylmm r nmA n <g> B m , where {A n } and {B m } are bases 
of Hermitian operators and r nm are the matrix elements of the correlation matrix R, the 
quantum discord D(A : B) is non-zero if the rank of the correlation matrix R is greater than 
the dimension of system A (similarly, D(B:A) is non-zero if rank(R) > dim(B)). That is 
to say that in order to show that a system has non-zero discord the number of measurements 
needed is only as many as is needed to show that rank(R) > dim(A). 

In our experiment, we find the rank by measuring several columns of R and examining 
the singular values to determine a lower bound on the rank. If this lower bound is greater 
than the dimension of system A, we can conclude that the system contains non-zero discord. 
If the lower bound is less than or equal to the dimension of system A, then we measure 
an additional column and find the singular values again. We continue this procedure until 
either the rank of R is greater than the dimension of A, or full state tomography has been 
performed, in which case this procedure produces no conclusive statement regarding the 
presence or absence of quantum discord. 

Experimental System: We implement this algorithm using a nuclear magnetic res- 
onance quantum information processor in the liquid state. NMR computation at room 
temperature is highly mixed and thus, is an ideal system for implementing the DQC1 model 
of computation. The state of the system in NMR is written as pnmr = (1 — a)I/2 N + ap pps , 
where p pps is known as the pseudopure state with unit trace and N is the total number of 
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qubits. Note that p pps need not be pure and the name is used for historical reasons. Mea- 
surements of ppp S are compared to a known state that serves as a reference for a. At thermal 
equilibrium, a is the ground state bias as given by the Boltzmann distribution and is equal 
to h'-fB /2ksT, where B° is the value of the static magnetic field, 7 is the gyromagnetic 
ratio, and /cgT is the thermal energy. For our experiment, with carbon-13 nuclei at room 
temperature in a 16.4 T magnet, the polarization is 1.4 x 1CT 5 . For this polarization and 
perfect implementation of the unitary transformations, the numerically computed discord 
present at the completion of the algorithm is 5.4 x 1CT 11 . Analytical results [3] for the 
average discord after a DQC1 circuit for a unitary drawn uniformly by the Haar measure 
indicate a discord of approximately 7.1 x 1CT 11 . 

The existence of non-classical correlations in an NMR state, as measured by the quantum 
discord, does not depend on the value of the polarization, provided it is non-zero - this can 
be deduced as follows: It has been shown [12] that a bipartite state has zero discord D(A : B), 
if and only if there exists a projective measurement {E^} on system A such that 



Noting that equation (|2]) is true for pmm if an d only if it is true for p pps , indicates that if 
non-zero discord is detected for p pps then the full NMR state also contains non-zero discord. 

We experimentally test for the existence of quantum discord in a liquid state NMR 
implementation of a four qubit DQC1 algorithm. Our qubits are the four carbon nuclei in 
the carbon-13 labelled molecule trans-crotonic acid (see Figure [2]). The qubits are a bulk 
ensemble of approximately 10 20 identical spin-1/2 nuclei that are manipulated in parallel and 
exhibit a two level energy splitting in the presence of high magnetic fields. Measurements 
of (a x ) and (a y ) are performed using quadrature detection of the free induction decay. 

The experiment is implemented in a Bruker Avance 700 MHz spectrometer where the hy- 
drogen nuclei are always decoupled using the WALTZ-16 and bi-WALTZ-16 composite pulse 
sequences [T3]. Radio frequency (r.f.) pulses in the plane perpendicular to the static magnetic 
field are numerically generated using the GRAPE algorithm [TH H5] , which starts from a ran- 
dom initial guess and is iteratively improved through a gradient ascent search. These pulses 
are optimized to start and finish with zero amplitude, produce a fidelity \to(Ug oal U S i m )\ 2 / d 2 , 
where d is the dimension of the Hilbert space of U goa i, of at least 0.998, and be robust 
against r.f. inhomogeneities (±3%). Before running the experiment, the pulses are adjusted 
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FIG. 2. (Color online) The molecule trans-crotonic acid and a table with the parameters of the 
Hamiltonian given in Hz. The shaded diagonal elements represent the chemical shifts with the 
Hamiltonian £j7TCo>j<7*. The remaining elements indicate the scalar coupling constants Jjj with the 



for non-linearities in pulse generation and transmission by placing a pickup coil inside the 
bore of the magnet and running a feedback loop, iteratively making small changes to the 
pulse shape to ensure as close to perfect transmission as possible. After the feedback loop, 
pulses are found to have a simulated fidelity of approximately 0.99 with respect to the ideal 
unitary. 

In this experiment, the bipartite split is the natural choice of the first, slightly po- 
larized qubit, and the remaining three. Thus, our state can be represented as pnmr = 
(1 — a)I m /2 A + o.i^2, nm r nrn A n <g> B m ), where only p pps has been written in terms of Her- 
mitian observables. We choose the observables to be {A n } = {I , X,Y, Z} and {B m } = 
{III,IIX,IIY,IIZ,IXX,...}, where X = a x (similarly for Y and Z) and {B m } are all 
possible combinations of three qubit Pauli operators. The size of our correlation matrix R is 
then 4 x 64 with a maximum rank of 4, and non-zero quantum discord D(A : B) is witnessed 
when rank(R) > dim(A) = 2. The unitary transformation used in this instance of the 
DQC1 model is U = diag(a, a, b, 1, a, b, 1, 1), where a = -(e^/ 5 ) 4 and b = ( e - i3 ^ 5 ) 8 , and 
is an important transformation in the approximation of the Jones polynomial for a class of 
knots whose braid representations have four strands. The problem of evaluating the Jones 
polynomial is not efficient using current classical methods, but has recently been shown to 
completely encapsulate the power of the DQC1 model P2] and has been experimentally 
implemented to successfully distinguish a specific class of knots p2] . 

Results: We tested for the existence of quantum discord at both the beginning and 
final stages of the DQC1 algorithm. To determine the rank, we choose to first measure four 
columns of the correlation matrix, denoting it Rtmnc, the truncated correlation matrix, and 
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test to see if the rank is greater than two. If not, we will measure additional columns until 
we find a rank greater than two or reach full tomography. To estimate the spread of the 
singular values as a result of uncertainties in the experiment, we perform a Monte Carlo 
sampling to determine which singular values can be reliably distinguished from zero. Shown 
in Figure [3] are the distributions of singular values for the (a) initial and (b) final states. The 
sampling results are binned and normalized by the total number of samples and bin width 
to produce a distribution of the relative occurrence. The cumulative of each distribution is 
included to guide the reader in estimating the integration of portions of the distributions. 

In the initial state, shown in Figure |3^, there is one clear non-zero singular value and 
three others that are very close to zero. Because the rank was not determined to be greater 
than 2, all columns of the correlation matrix were measured. The histogram is generated by 
sampling from 1000 random combinations of four columns, always including the first column 
and sampling from their errors ten times each, giving a total of 10, 000 samples. 

To measure the final state we started by measuring the columns corresponding to B m = 
III,IZI,IIZ,IZZ; recall that the four rows refer to A n = {I , X,Y, Z}. The truncated 
correlation matrix is shown in equation ([3]), and for example, the (2, 3) element corresponds 
to the observable XII Z and is equal to —0.13. Each matrix element is calculated using 
linear inversion of variables fitted directly from the NMR spectrum. In order to obtain these 
results, four instances of the experiment were required, each with a different readout pulse 
and observing a different spin. The uncertainties reported in equation ^ are propagated 
through linear inversion from spectral peak-fitting, and correspond to 68.2% confidence level 
in the results of the fitting process. 

/ 1 -0.01 0.00 -0.01 \ 
0.10 -0.34 -0.13 0.25 
0.17 0.38 0.04 0.26 
\0.01 0.08 -0.01 0.02 J 

( 0.007 0.01 0.007 \ 
0.05 0.05 0.05 0.05 
0.05 0.05 0.05 0.05 
\0.04 0.007 0.007 0.007 / 
Note that the IIII term has no error associated with it as it is assumed to be exactly 1 
since tr(p) = 1. All values in the correlation matrix are measured as a fraction of a by 
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comparison with a known reference state. Upon inspection of the results matrix in equation 
(|3j) it can be seen that there are likely at least three linearly independent columns (columns 
1, 2, and 4) due to the relative signs of the largest elements. This is confirmed by examining 
the singular values of Rtmnc, the distribution of which is shown in Figure [3}d. Because three 
non-zero singular values were found immediately no additional columns of R were measured. 
Therefore, non-zero discord is reliably detected in the final state of the DQC1 algorithm, 
even though simulations indicate the amount of discord present is 5.4 x 10 -11 . 

Conclusion: Since the polarization in the DQC1 algorithm does not affect the presence of 
discord in the full quantum state, we successfully detected a very small amount of quantum 
discord in a bulk ensemble, highly mixed, liquid state NMR implementation of the DQC1 
algorithm. There was no conclusive evidence for the presence of discord at the beginning of 
the computation, consistent with simulations indicating discord is generated in the DQC1 
algorithm. Whether or not these correlations are of use in the computation is a question of 
great interest and will shed more light on our understanding of quantum systems. 

G.P. would like to thank M. Ditty for his technical expertise with the spectrometer. This 
work was funded by NSERC, Quantum Works, and CIFAR. 

Upon completion of this manuscript, we became aware of a concurrent independent 
study [18] that corroborates our conclusions, but uses a different method. 
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a.) Initial State 
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Singular Values 

FIG. 3. Shown are the distributions of the singular values (and their cumulative distributions) 

computed for the experimentally determined correlation matrices of the (a) initial and (b) final 

states of a DQC1 algorithm implemented in NMR. These distributions are created by sampling 

from a normal distribution of the errors on each matrix element and calculating the singular values 

of the sampled matrix. There are 10,000 samples in each plot and the histogram bin size is 0.005. 

The cumulative of each distribution is included to guide the reader in estimating the integral of 

portions of the distributions. We deduce that the correlation matrix of the initial state has a rank 

of 1 and the correlation matrix of the final state has a rank of at least 3, which is to say that we did 
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not detect quantum discord in the initial state, but did witness the presence of quantum discord 



